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1. Introduction and Main Theorem 

In his talk at JAMI 2009 on March 24, Kurokawa presented a rather mysterious 
looking identidity of elementary number theory: 
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Example 1.1. n = 12 = 2 2 • 3: 

— (1 + 2 + 3 + 4+1+6 + 1 + 4 + 3 + 2 + 1 + 12) = ^ = y 
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In fact, this is a special case (the case r = 1) of identities obtained by [KO| : 



Theorem of Kurokawa-Ochiai [KOj . For n, r £ N, 

i ™ r r / i 

(2) — gcd(n,fci---fc r )=JJ E ordp( „)^ M-- 

fci,...,fe r =l p|n U=0 ^ ^ 

Here, m H r is the repeated combination: 

m H r := $ {degree r homogeneous monomials in m variables} 
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Example 1.2. n = 6 = 2 • 3, r = 2: 
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--(1 + 2 + 3 + 2 + 1 + 6 + 2 + 2 + 6 + 2 + 2 + 6 + 3 + 6 + 3 + 6 + 3 + 6 
6 Z 

+ 2 + 2 + 6 + 2 + 2 + 6 + 1 + 2 + 3 + 2 + 1+ 6 + 6 + 6 + 6 + 6 + 6 + 6) 
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[KO obtained © by studying some multivariable zeta function of Igusa type, 
and Kurokawa said he is not aware of any elementary proof even for (JT|) . 

Now the purpose of this paper is to give a purely elementary proof of general- 
izations of Kurokawa-Ochiai identities (JSJ) from the view point of elementary prob- 
ability theory. Fixing n £ N, we would like to understand the random variable: 



(3) 



X : n := N r N 
(ki, k 2 , ■ ■ . , k r ) h-> gcd(n, kik 2 ■ ■ ■ k r ), 

Although fl = N r is an inhnite set, we would like to regard it being equipped 
with the "homogeneous measure" . For this purpose, we observe: 

ki = k[ mod n (i = 1, 2, . . . , r) =>■ gcd(n, kik 2 ■ ■ ■ k r ) — gcd(n, fcifc 2 ' ' 

Thus, instead of ([3]), we may equally consider the following random variable: 

X : Q := {1,2, . . . , n} r — > N 

(ki,k 2 , . . . , fc r ) i-> gcd(n, kik 2 ■ ■ ■ k r ), 



K) 



(4) 



where O = {1,2, . . . ,n} r is equipped with the homogeneous measure. Then the 
Kurokawa-Ochiai identity ([2]) is nothing but a convenient formula to evaluate the 
average E[X] of the random variable (j4|). Now, from the view point of elementary 
probability theory, it is very natural to seek for similar convenient formulae for the 
variance V[X] = E[X 2 } — i?[X] 2 and even "higher" invariants. 

Our Main Theorem offers such formulae for the continuous version E[X W ] (w £ 
C) . Their special cases w £ N are nothing but the moments of the randowm variable 
X, and the simplest case w = 1 is nothing but the Kurokawa-Ochiai identity ([2|). 

Main Theorem. For n, r e N, w G C, 



22 gcd(n, k! ■ ■ ■ k r ) 1 



Ylp\n Sz=0 ord p (n)Hl (l ~ ~ 

-I lp] n 
1 Ip| n 

Corollary 1.3. For n G N and w 6 C satisfying w = 1 or |p 
(5) '■ 1 



1-P" 



l_ p w- 



(1 



i/iu / 1 
ifw = l 



(^r) r +^ (n)( - 1) (i-P" 1 rESH(.)fli {(i-P _1 ) i_r - (i -p"'- 1 ) 1 -'"}] « # i 

EL ord p (n)^i (l - ~) iflV = l 



1| > 1, 



lim 



gcd(n, fei ■ • ■ fc r ) 



fci, 



Heuristically speaking, ([5]) implies the random variable ([3]) "converges" to the 
constat random variable concentrated at n. 

When n is fixed, the number of the terms needed to be evaluaed in 

Sfci k =1 § cc K n ? ki ■ ■ ■ k r ) w grow exponentially with respect to r, but Main 
Theorem allows us to evaluate it with lear growth with respect to r. 

In our purely elementary proof for Main Theorem, the original identity {T]) turns 
out to be a matter of triviality. In the course of proving the general case, we shall 
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show some identitites involving the repeated combinations (Xcmma ^.ip . which may 
of of independent interest. 

We note another kind of generalizations of |T]) is given in [DKKj . which is once 
again proved using some zeta function of Igusa type. Even for this and some 
generalization, we can offer an elementary proof |M1| . It is possible to define and 
study some multivarible deformed zeta function of Fi-scheme of Igusa-type [M2j , 
which generalizes both |DKK] and |KOj . 

The basic idea of the results in this paper were obtained during the author's 
stay at JAMI2009, Johns Hopkins University, in March 2009, and presented at 
NCGOA2009, Vanderbilt University, in May 2009. The author would like to express 
his gratitude to Katia Concani, Alain Connes, Jack Morava, Takashi Ono, Steve 
Wilson, and Guoliang Yu for their hospitalities. The author also would like to 
express his gratitude to Nobushige Kurokawa for his work and encouragement. 

2. Proof of Main Theorem 

Whereas [KO] used some multivariable Igusa type zeta funciions for groups, we 
exploit some ring structure: For a finite ring with n elements R = {ki G R | 1 < 
i < n} and r G N, set 



ZrH :- 



1 



\R\ 



Oi,...,fc r )efl r 
We wish to understand this, because 



E \R/{k 1 ---k r )R\ w (iugC) 



1 n 

fel,...,fe r =l 

Of course, we have an elementary probability theoretical interpretation: For the 
random variable 

X : n := R r -> N 
(fci,fc 2 ,...,/c„) h-> \R/(k v ■ -k r )R\, 

where £1 = R r is equipped with the homogeneous measure, 

E[X ' W] =\W £ \R/(kvk r )R\ w = Z r R (w) 

' ' (k!,...M r )eR r 
Elementary properties of Z r R (w): 

• If there is a finite ring decomposition R = \\ i Ri, then 

zu™) = Y[ z rA w ) 



. then 



Set N R (f) := | {(ki, ...,k r )€R r \ \R/(h ■ ■ ■ k r )R\ = /} 

\R\ 

ZrH := ^ E \ R I^ ■ ■ ■ = ^pE N nWf 



{k 1 ,...,k r )GR r 



/=1 
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So, if n — Y[p\nP m ' dp ^ > then by the Chinese Remainder Theorem, 



1 " 

- gcd(nM---k r ) w =Z r z/nZ (w) = l[Z r z 



Z/p ord p(">Z 



(w) 



(0) 



fci,...,fc r =i 



rdp(r>) 



II „rord p (n) E ^Z/p«z(/) 



. p-/ uidp(n) 
p|n /=1 



Thus, we need to compute 

%. z (/) = | {(fci,...,^) € (Z/ P e Zy\\(Z/p e Z)/(k 1 ---k r )(Z/p e Z)\ = /} 

= {(fci, . . . , fcr) € {1, 2, . . . ,p e } r \ gcd(p e , h ■ ■ ■ kr) = /} 

for 1 < / < p ord p("). 

To get some feeling, we first play with the simplest case r = 1, Then, it is easy 
to evaluate: 



p£—d _ pC-d—l 



(7) 



^/ p e Z (/) = { 1 





if / = p d (0 < d < e) 
if /=p e 
if otherwise 



which follows from 



# 



{1,2,. ..,p e } 



ord„— c£ 



■ p 

ord„ 



/=1 



Proof for r = 1: 

fc=i 

TT 1 

X J. ^ordj 
P I '< 

=n 



pjn 
ordp (n) — 1 



p\n 



ordp (n) — 1 

(i - p- 1 ) (p w ^) d \ + p° rd * (n)(, °- l) 



d=0 



n p |Jl+(l-|)ord p (r ' 



if it; = 1 

if to ^ 1 
if w> = 1 



□ 

So, the case r = 1 is a matter of triviality, and not much extra effort in consid- 
ering general w G C. However, the proof for general r G N is more complicated. 



n p |n(l+(l-|)ord P ' 
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Observe for < d < e, 



{(P e -P e 1 )-l + (p e 1 -P e \-(p-l)x e } r 

= (p e - p e ~ 1 y x The coefficient of x d in ^ 1 + f + • ■ • + ( 



V ) VP 



(8) 

= (p e — p e ) r x The coefficient of £ in < 1 



e-l , 

= (p e - p e_1 ) r x The coefficient of x d in ^ r #, I - 



= (P e -P e ^''rHd 



To evaluate N^ pl , z (j) e ) and for further computation, we prove the following iden- 
tities involving the repeated combinations, which may be of independent interest: 

Lemma 2.1. For e,r £ N, set a degree e — 1 polynomial of x by 

e-l 



(9) f e {x) :="£ r H k x k 

k=0 



Then, 

(10a) /T(l) = r+l-^e— 1 = e+r-lCe-l = e+r-l^r = e^r 

(10b) (1 - xYfM + y e r;(l -x) = {l- x) r (r e (y) - y e /;(l)) + y e f r+1 {l - x) 
(10c) (l-xyr e (x)+x e fe(l-x) = l 

Proof. (jlOap is an easy consequence of the Taylor expansions at x = 0: 

e-l e-l e-l 

/J(l) = rH k = The coefficient of x e ~ x in ( ^ r i/ fe 2; fe ) ( ^ x fe ) 

fc=0 fc=0 fe=0 

= The coefficient of x e ~ l in ((1 - x)~ r {l - x)' 1 = (1 - x)" (r+1) ) 
— r -\-\H e —\ — e j r , r —\C e — i — e j tr —\C r - — cH r 

Now (|10b[) is an immediate consequence of f)10a|) : 

(i - x yrM + y e m - x) = (i - s) r (/ e r (i/) - z/7;(i)) + z/ e (iV e (i - x) + - x y) 
®m (i _ x y ( /:(2/) _ y7:( i)) + 2/ e ( (E eH k (i xA + Mi xy) 

(i - x) r {rM - y7e r (i)) + y e Eeff fe (i - x) fc i (i - a;) r (/!(?/) - y e /e r (i)) + 2/7 r e +i(i - *) 
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To prove (|10cp . consider the Taylor expansions of degree e of (1 — x) r at x = 0, 
and degree r of x~ e = (1 — (1 — x)) e at x — 1, respectively. 

e-l 

(11a) f e {x) = rH k x k = (1 - x)- r + x e ■ 0(l;x - 0) 



fc=0 
r-1 



(lib) #(1 - X) = eH k (l - X) fc = X~ e + (1 - x) r • 0(1;. T -> 1) 

k=0 

We now set the left hand side polynomial of x in (J9J) to be 
(12) l(x):=(l-x) r f r e (x)+x e f?(l-x) 
Then, applying (|lla|) and ()llb[) respectively to (fT2j) . we see 

l(x) = {1 - xf ({1 - x)- r +x e -0(1; i +x e /;-(l-^) 

(is) ; 7 x 

= 1 + x e ■ ((1 - z) r O(l; a; -> 0) + / r e (l - a:) J 

i(ar) = (1 - x) r f r e {x) + x e (x- e + (1 - .T) r • 0(1; x 1)) 

(14) ' V , 7 

= 1 + (1 - x) r • (/J^) + .T e ■ 0(1; x -> 1)J 

(fT5| and (Q3J) respectively imply i(x) — 1 is divisible by x e and (1 — x) r . Thus, 
l(x) — 1, a polynomial of x of degree at most e + r — 1, is divisible by x e (l — x) r . 
Of course, this implies l(x) = 1, and thus ()10cp has been proven. 

□ 

Now we may complete our computation of N^ pez (f): 

Proposition 2.2. For t/ie cyclic ring Z/p e Z p a prime and e G N, 
(15) 

f( p e_ p e-l ) r rifd ^ d if f = P d (0 < d < e) 

Nl /rZ {f) = p e(r-l) f e (l - l) = p*-D ^ ^ (l - l)' z/ / = ^ 

[ i/ otherwise 

Proof of Provosition \2.2l It is obvious that N^ peZ (f) = if / is not of the form 
p d {0 < d < e). The case of < d < e is already treated in ([5]). 
Finally, the case of d = e can be taken care of as follows: 



d=0 d=0 



(pt - (p e - P e - L ) r f r e ( - j = p er 1 1 - ( i - - j /; 



^ r ( - ) /; [i--)= P e(r - 1} /; [i--)= P e(r 



v j \ v i \ p 



X>(i-;)' 

(=0 v ^ 7 



□ 
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Proof of Main Theorem. 



k 1 ,...,k r = l p\n V f = i 

W TT 1 

1 1 jf ordp(n) 



ordp (n) — 1 



=n 



ordp (n) — 1 



n 



/t —l\r pr f w~ \\ , ( w— l\ord p (n) „ordp(n) /l — 1\ 

(l-p ) /o,d p („)(p ) + (p ) /r * v ) 



flfibl 



1 lp| n 
1 lp| n 

rip|,i 

1 Ipln 



(!-P ) /ord p (n)(P ) + (P j fr ) 

(l-P ) (/ord p (n)(P )-{P ) /ord^njl 1 ) J + [P ) fr + 1 (1 ~ P ) 

if 10 # 1 



if u> 
if to 



/^"'(I-P- 1 ) 



if 1/7 = 1 



ffl f IU [(^)" +P ordp( " ,(m - 1) ESo^wfl^i-Pl' - (^) r (i - p- 1 )'}] if « # i 

|n p |„ E!=0 ord p (n)Hl (l - ^) if W = 1 

'ru [(r^r) r + P ord " ( " )(i "- i) (i -p-'rEiToord.w^ {(i -p- 1 )'^ - (i -p 1 - 1 ) 1 -}} nw^ 

Il p \n E[=0 ord p (n)Hl (l - | ) if W = 1 



□ 



Proof of Corollar il.3\ When it) = 1, 



lim TT 

-»+nn J- i. 



r — >+oo 



p\n 11=0 



y^ord p (n)-^ I 1 



n 



+00 , 1 

X^ord p («)-ff; I 1 - - 



"^'nfi-M 



- ordp(n) 



ordp (n) 
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When \p w - 1| > 1, since |1 -p^\ < 1< |1 -p w \, 



lim rr 

n 



i- P - 



. ordr, (n) (iv — 1) \ ^ rr I /i 



lim 

r — > + oo 



1=0 



l\l ( 1 ~P 



n 



i -p™- 1 

rr£r) +^ ( " )K1, E^^.(i-^)' 

r-1 

p p{ A ; hm (1-p ) > ord p („)Hz(l -p 

z=o 

J^J j^Q _^p° rd p(^)(^- 1 )p° rd p(^) 

p| n 

_ p o rdp (n)( W -l) rj (1 _ p -l)rg (n)ffi (l_p— 1)' 

Z=0 

r-1 

«d, ( n)w_ ori,(„)(»-i) lim (i-p-YV^^ji/^i-p- 1 ; 



1 -p ) 



This is shown to be equal to n w , if the following is shown: 

r-1 

lim (1 - V ord p (n)^(l - p W - r ) l - r = 



However, this follows from 

r-1 



lim 

r — >+oo 



{i-p- x ) r Y,°^(n)m-p w - l ) l - r 



1=0 



= lim 

r — >+oc 



E(=0 ord p (n)Hl( 1 _ pW -i ) 



(I-P" 1 ) 



n-i 



< lim 

r — >+oo 







1 v-' 




(l-p-l)-l 


r 



r — t'+oo 



(i-p- 1 ) 



n-i 



® ,. Jr (1) {HP ,. ord p (n)+r-l^ord„(n) . 

hm - = hm — f 77 — = 0, 



11-1 



1\-1 



where the last equality follows since ord p (n)+r-iCord p (n) is a degree ord p (n) poly- 
nomial of r and |(1 — > 1- Now the claim has been proven. □ 
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